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(57) Abstract 

An ultrasonic measurement system performs a signal path measurement by directing ultrasonic signals through a gaseous material 
in a conduit, and processing the detected signals to determine sound speed and to derive the average molecular weight of an unknown 
hydrocarbon mixture present in the material. The processor includes a plurality of stored tables of critical constants of hydrocarbon mixtures 
as a function of the average molecular weight of the mixture, and is configured to iteratively set a hypothetical molecular weight, determine 
the corresponding critical properties, and compute a predicted sound speed. If the two speeds differ, a new weight is set and the procedure 
is repeated until the predicted sound speed matches the measured speed, indicating that the current estimate is the correct average molecular 
weight Once the processor has determined the critical constants from its stored tables, it applies the virial equation of state and mixing rules 
to determine the predicted sound speed for each hypothetical molecular weight of an unknown mixture of hydrocarbons present, together 
with one or more known inorganic components, in the fluid material. The processor may be configured for user input of known quantities 
of one or more gases such as nitrogen, hydrogen, sulfur dioxide, carbon dioxide or other inorganic or non-hydrocarbon gases present in the 
conduit. Alternatively, all or some of this data may be provided in an automated manner from suitable signals, settings or measurements 
from upstream conduits, valves of measurement and control instrumentation. 
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ULTRASONIC MEASUREMENT SYSTEM 
WITH MOLECULAR WEIGHT DETERMINATION 

Technical Field 

5 The present invention relates to ultrasonic measurement systems, for example 

systems such as are used to determine properties of a fluid flowing within a conduit. 
These systems operate by creating an ultrasonic signal burst, and transmitting the burst 
into the fluid and receiving the signal after it has traveled along a path through the fluid. 
Two or more related paths may be used (for example upstream and downstream paths, 

10 orthogonal paths, or clockwise and counterclockwise paths) to develop signals indicative 
of sound speed, or of directional differences in sound speed that depend on flow 
components. Generally the property of the fluid is determined as a function of the transit 
times of these signals. 



15 



20 



25 



Often it is desired to obtain from the signals useful information such as the mass 
flow rate, or the caloric value of the material flowing in the conduit. In that case further 
information is necessary and this often requires the system to operate in controlled 
environments where the user enters extrinsic data via a keyboard to specify the known 
components of material in the flow so that the signal processor can perform its 
calculations with the corresponding values of physical parameters for those materials. 
However, such operation is only possible in a limited number of well controlled 
circumstances. In some plants, and especially in situations where flue gases are present, 
the composition of the gas stream will vary over time. It may also be subject to certain 
fast fluctuations, as occurs in flare gases or certain hydrocarbon processing 
environments. When the gases present under these measurement conditions are used in 
bulk processes involving hydrocarbons, the monetary value of material flowing in the 
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conduit can be substantial, and it may be necessary to know its volume and composition 
with relatively high accuracy lor, processing, cpntrol or environmental reasons. 

It would therefore be desirable to extend the u?^ of 
ultrasonic measurement systems to include flows with varying but unknown composition 
of hydrocarbon components. - • ,- ■ « , 

It would also be desirable toprovide an ultrasonic measurement system which 
determines the molecular weight of a component of a fluid flow. 

Summary of the Invention 



These and other useful ends are achieved in accordance with a basic embodiment 
of the invention wherein an ultrasonic transducer propagates a signal through a fluid and 

15 determines the average molecular weight of an unknown mixture of hydrocarbons present 
in the fluid. A processor receives sound speed, pressure and temperature information, 
and performs an iterative procedure to match a hypothetical composition or average 
molecular weight to the observed sound speed. In accordance with a preferred aspect of 
the present invention, the microprocessor stores a number of tables in which the critical 

20 pressure, compressibility, acentric factor and heat capacity of hydrocarbon gas have been 
tabulated as functions of molecular weight. The microprocessor executes a programmed 
set of calculations for determining the sound speed of a mixture of gases by using the 
virial equations, determining the virial coefficients from these four properties of the gas. 

25 The processor includes a plurality of stored tables of critical constants of 

hydrocarbon mixtures as a function of the average molecular weight of the mixture, and 
is configured to iteratively set a hypothetical molecular weight, determine the 
corresponding critical properties, and compute a predicted sound speed. If the two 
speeds differ, a new weight is set and the procedure is repeated until the predicted sound 
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speed matches the measured speed, indicating that the current estimate is me correct 
average molecular Wight! Once the processor has determined the critical constants from 
its stored tables, it applies the virial equation of state and mixing rules to determine the 
predicted sound ^ s^^#each '^tnetical molecular weight of an unknown mixture of 
hydrocarbons ^ present; together Vith one or more known inorganic components, in the 
fluid material. The processor may be configured for user input of known quantises of 
one or more gases such as nitrogen, hydrogen, sulfur dioxide, carbon dioxide or other 
inorganic or non-hydrocarbon gases presentin the conduit. Alternatively, all or some of 
this data may be provided in an automated manner from suitable signals, settings or 
measurements from upstream conduits, valves or measurement and control 
instrumentation. The hydrocarbon analysis may be used to detect leaky valves or 
upstream process occurrences, or to evaluate gases in feed loops or outflow condmts. 

Brief Description ot' the Drawings ■ 

' These and other features of the invention will be u*ierstood from the description 

herein, taken together with the Figures wherein 

Figure 1 mustra.es a reiationship between molecnlar weight of hydrocarbons and 

acentric factor; 

Figure 2 illustrates a relationship between mofceuiar weigh, of hydrocarbons and 

ideal gas heat capacity at 220°K; 

Figure 3 illustrates a relationship between molecnlar weigh, of hydrocarbons and 

ideal gas heat capacity at 298 °K; 

Figure 4 illus.ra.es a relationship berweeu motecular weigh, of hydrocarbons and 

ideal gas heat capacity at 560°K; '• • • 

Fignre 5 illustrates a relationship between moleraiar weight of hydrocarbons and 

critical temperature; " : ' '' '"' ' ' : ' 

Figure 6 illustrates a relationship between n^lecular weight of h^^^^ 

critical volume; 
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Figure 7 illustrates a relationship between molecular, weight of hydrocarbons and 

critical compressibility ; and . . /i,- ^ha h 

Figure 8 compares sound speed, in ethylene as calculated in accordance with the 

invention, arid as measured; . . . :j : ; .. , 

5 Figure 9A and 9B compare calculated and measured sound speed for two natural 

gases as a function of pressure; , ; , , , K . .- , 

Figure 10 compares errors using the molecular weight determination of the 
' invention with that of ajcndwn protocol for a range of . gases; 
Figure 11 illustrates operation^ of jthe invention; and 
10 Figure 12 compares predictive calculations to empirical properties of six gases. 

Detailed Description 

The present invention relies upon Applicant's discovery that for a broad class of 
15 hydrocarbons , each of the critical phy sical properties of the gas which are necessary for 
performing the virial calculations is itself a linear or. very well behaved function of the 
average molecular weight of the hydrocarbon mixture. Applicant realized that, by 
compiling a calibration table for a number of different mixtures spanning the expected 
range of molecular weights, one can accurately determine all four of the physical 
20 constants required for virial calculations from this single value, the average molecular 

weight. In particular, applicant has found that for common hydrocarbon gases including 
aromatic, double bonded, and single chain hydrocarbons, the relationship is substantially 
independent of the particular amounts and types of each chemical structure in the 
mixture, and depends simply on the average molecular weight. Thus to compile and 
25 store reference tables for converting a molecular weight into .critical temperature, critical 

pressure or volume, critical compressibility , and acentric factor, it is sufficient to perform 
these measurements for a large number, e.g. , 50-150, of different mixtures of 
miscellaneous hydrocarbon gases that span the expected range of average molecular 
weights to be found in the field. Once this is done, gas properties may then be modeled 
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in terms of a single variable; the average molecular weight. Apparatus according to 
applicant's invention includes a processor which, by simply changing this, single variable 
and comparing calculate^^iiiionies with the actual measured outcome for a simple 
property, such as sound speed, performs an iterative procure to determine the actual 
average molecular weight of mdfcyarocWbda fiuxture. • This is a measured of 
importance in many petrochemical processes, and in further embodiments allows one to 
directly calculate caloric' value of the gas or locate sources of the component 
hydrocarbons entering the upstream flow, the process also yields other properties of the 
gas, such as the critical compressibility, useful for process control or for deriving a mass 
flow or other measurement form the sensed parameters. 

Figures 1-7 show the relationship between average molecular weight of a mixture 
of hydrocarbon gases with other background gases, and the gas properties of the mixture. 

Skipping ahead bfiefly to Figure 1*. the iterative procedure proceeds as follows. 
The system is configured for a user to enter the percentage and types of inorganic gases 
which are present in the mixture such as hydrogen, nitrogen, sulfur dioxide, carbon 
dioxide and, as needed, any other commonly present inorganic gas. This information 
may be user-entered, if known, or provided from various sensors and controls, such as 
pressure and flow transducers or valve setting or control signals from the upstream lines, 
to thediscussion beioW; it wUl be assumed that the relevant physical constants for these 
background gases, all of which have' been extensively studied and^tabulated in the 
literature, are stored in the processor, and that if a non-standard background gas is 
involved, the necessary constants have been determined and stored. The processor then 
proceeds to estimate the molecular weight of the hydrocarbon component of the fluid 
stream, and perforin a sequence of estimates and calculations to result in a predicted 
sound speed, whicfi is checked against the measured sound speed.. If different, it sets the 
next molecular weight estimate, and checks its effect on predicted speed, until these 
converge to the average molecular weight of the hydrogen-carbon components present in 
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the gas. '• v. . ■*.•='>.•,: ^..\ 

This is done by first setting the molecular weight Jow, e.g., equal to that of 
methane. Using its stored tables the processor determines four values (critical 
compressibility, acentric factor, and so forth) conespopding to that molecular weight, 
and using the measured temperature, pressure and sound speed applies the virial 
coefficient method and mixing roles to determine the predicted sound speed for a 
hydrocarbon gas with those constants (determined from the set molecular weight) mixed 
with the amounts of the specified morgan^ gases known to be present ^ 
may, for example, be heuristically known or measured. 

The sound speed predicted by the virial equations and mixing rules is then 
compared to the measured sound speed, and if the results differ, the processor again sets 
a molecular weight, for example, at the other extreme of the range, e.g., sets mw=130, 
and again goes through the process of looking up the critical properties corresponding to 
the tentatively set molecular weight, applying the virial equations to determine an 
estimated sound speed of the set mixture, and comparing the estimated sound speed with 
the sound speed actually measured by the transducers. Once the molecular weight has 
been determined, the output may be used to set desired flow or recirculation based on the 
caloric value (when used as a combustion feed), or may be used by a decision module or 
control device to identify the upstream source of one or more components in the flow. 
This is especially useful to identify aberrations such as stuck or leaky valves giving rise 
to the presence of an unexpected outflow component. 

The microprocessor proceeds in this way, interpolating or extrapolating from 
successive molecular weight estimates until the set molecular weight yields a theoretical 
sound speed which matches the observed sound speed. This value is then taken as the 
average molecular weight of the hydrocarbons present in the flowing material. A more 
detailed discussion of the relevant virial equations and mixing rules follows, illustrating 
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representative calculations to demonstrate the accuracy achieved by the foregoing . 
method. Because the method employs a fundamental relationship, namely the virial 
equation of state; and that equation and mixing rules have been w^U verified, the output 
of the processor imy approach me best accuracy possible -for this type of calculation on a 
complex and unknown mixture. < ■ . , • ; 

A representative set of processing steps are discussed below. Background 
discussions of the virial calculations and the mixing rules for calculating the effects of 
mixing several gases together may be found, in texts and common treaties, for example at 
chapter 3, section 11, chapter 4 of the Properties of Gases and Liquids, 3rd Ed., R. Reid 
et al. McGraw-HUl Book Company. Applicant has tabulated the acentric factor, heat 
capacity at 220°K, 298°K and 560°K temperatures, critical temperature, critical volume 
and critical compressibUity of hydrocarbon mixtures as a function of molecular weight 
and these tables appear as FIGURES 1-7, respectively, herein. The Figures further 
show the corresponding values of these properties for hydrogen, nitrogen, carbon dioxide 
and sulfur dioxide, by way of illustration for inorganic background gases. 
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Ekjuaticmsfor <0*3 VirialMetLLod 

. _ , / T o Hirschfeider, C. F. Curfciss, and 

The virial equation for sound speed is CJ -V-/" John Wiley and Sons 
R. B- BirdKoieca'tor TVieory of. Gases and ; ^7"^~ * i ^ ^, ^ 
1964.P-232) 



1. 



10 



15 



. . , , i,. c p -a in the *ras constant, T is. the _ ■ 

where y 0 is ^ idead gas c~ ' ' ' 

M is the molecular U>4 V is the tote M B 

* thLoodvhial coefficient, * units of **™£ZtL d!U 
Thesccoaoviriarc^t^ 

temperature, critical pressure and the m aaaified to fit 

F.Curl, JAOS79.2369^19SU «» *™SSS£^S: ^ Sherwood, 2fte 
noupotar molecules G. Reid, J. M. P«"S^^ 7 p 63 8S) . 
Properties of Gases and iiguKis . McGraw mu .w 

;1 + <-.13?5 + 331;w )• ^) 

* <-*!» - *** )• (t) 3 + <-- O0Ofi0 '" •° 0$ ' M V 8 
The molar volume V is given by 

-rr_R-T 



R-T, 



[ 



c I -x. /- 1385+ -331?<» )• 
( .1445 + .0537- co ) + -33— | + < -M" 



20 



P 



For convenience 
belov/ 



, the derivatives of B ^ <he above equation are listed 



R.T r T 7 T A 

g 3 - c - 33._£ + <^77--662«<a.)-l-=-J •» 



25 



+ ( .0363 + 1.269 



+ (.004856 + .064 



"I 2 * 



-.66«— + (--831 4- li966 <a ?>•« *" 
+ (-.1452- S^O-d) + (^3704- 



30 
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The virial coefficient for mixtures is given by. the expression 



where B £ j is the ^ ^cffidfentfor the binary miKture^f ^orients , .. 
i and j, and ,B it s is just, the virial coefficient of the pure component i. : 
X- and X- are mole fractions. ' < '• - " - 

Virial . W t*^TkS^M£S be 
constants. But as these are not generally toovm "ir ^ 

estimated using mjxing rules: 



C; S J C; * C; ^l.j 



4. 



Note that for single component data 



p =- — Y c p 



f : ~*>"i 9 flnd3, the pure component virial 

For example, for a ternary mixture 1, A J^£^ ^ £. critical data, using 
15 coefficients B j, and B 3 may be calculated ^ estimated ftom 
equation 2. The critical data for the 3 m^^vSal ^ 

Rations 4; and these date can he ^^^f s S^h equation S 
Sclents Bj2 , Ba3 ; and ^.f^^^^Fte this example 
to calculate the actual virial coefficient for the mixture, uw 
equation 8 becomes 

where and X 3 are the mole fractions of the three component,. 

E _ L£= rt r i for Determinme Mototo^lgM 

The goal is to obtain molecular ^^™g^&5?^ , 
25 gasef. The sound speed, ^?^^» ^ e ^^^ are presumed 

25 concentrations of any aoa-hydrcc^ons^ ^^pute molecular weight 

to be fcnovmu This is not ^^^^^gations, heat capacities, and 

«mnrr the eauations in part A above. ooncenaro.ww 

critical late are missing for the hydrocarbons. 
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The following scheme is used to estimate the missing data. Co^llMM 
are foundbetween molecular weight, and heat capacity and critical data 
for the hydrocarbons. For a presumed value of hydrocarbon molecular 
weight, these correlations provide values for j and© ; and derivatives. Ah^ 
overall virial coefficient is calculated using equation 4, and a sound speed is 
calculated using equation 1. The calculated sound speed is compared to the 
observed sound speed, and the calculation is iterated unfertile value of 
molecular weight which gives the observed value of sound speed is found. 

Iteration is not an instructive approach for this Mathcad document 
Instead, the sound speed will be found using the vmal equation for the 
known case of 85% methane in nitrogen. Then the vinal method will be 
anolied to determine what molecular weight gives the same sound speed, 
for the unknown case of 85% unspecified hydrocarbonCs) in nitrogen. 

Let X H :=0^5 be the hydrocarbon concentration, X B := 0.85 be 

the methane concentration, and X A ^= 0.15 be the nitrogen concentration 

The unknown molecular weight will be found to be roH = QM59 kg/mole. 

The temperature is T := 298.15 kelvins, and the pressure is P :=20 

atmospheres. 

1. Coefficients 

The correlation for heat capacity is determined using a weighted quadratic 
fit. The results are given in terms . of absolute temperatare T ^.^ e 
molecular weight nf H (kg/mole) of the hydrocaAons. Gv is m umts of joules 

per mole, with Rg as =8.S1441joules*mole"^^Sr 1 . 



R ^-.= 8:31441 



The hydrocarbon heat capacity coefficients are 



* 00 


:= 1.8414 


a 0i 


a 10 


:=-30348 


*il 


*20 


:= -682.06 


a 21 



f 3 E.Q2 I.S128- io" 6 
a ^1= -2.4852' 10" 5 
& 22:"= -5.7596- IO*" 3 



Ideal gas heat capacity C^B. for the hydrocarbons is given by 



c vH := R ^ 



(&oo +a or T+a o2'' ia ) - 

+ U 10 + a lt -T+ a n -T^ -mH ... 
+ (a 2Q-1- aarT+a^-T^-mH 7 



in units of Rgas 
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The other hydrocarbon data are the critical temperature T c (K), critical 
volume Y c 'lc^mpi^['c^iticai, Mmpressibaity and acentric factor <o. These 
are independent of . temperature! (^efficients "were ddtenmned fix»m ieas't 
square, fits. . •.. "., \ ~. X ,Y- ] . "_ r / [ .'..•.'..-} ... ,'. " .. 

c 0 := 41.053 0^-= 3080.8' ? c 2 608 1.2 

d 0 := .28566/ ; V" d [ :=' ; -?;i^443 r _ . ~ ; f ' : " . 

, c 0 :=3.ii4-io~' 2 :.cj .:= 16938 '•' l ' : ".'7V' ' 
The values are calculated from ^the -folio wing (mH in units of kg/mole) 

:=fa 0 +■ fa £-mH+ b 2 -niH 2 + b 3 -mH 3 
V ^ ;= c o +■ c aiH + c 2 - mH? 

r cH :=d o +d r mH 

Data for non-hydrocarbon gases are given directly. For example, for 
nitrogen (let nitrogen be gas A) 

m A := 0.028013 . := 126.2 v cA : = 89.645 

<a A :=0.04 \ . a 0A :=7.44 

a^^-3.24-10*" 3 ' a2A :=6.4vl0" 6 . a 3A :=-2.79- 10~ 9 

The coefficients for heat capacity are given as they are usually found, 
for C p in units of cal^mole"^^' 1 . For our uni^s 

C vA :=4.184- (a^-h a T+ a^T* + ^A** 5 ) - R gas 

using the conversion factor 4.184 joule/caL 

For purposes of demonstration only, let gas B be methane. 

m B := 0.016043 T <B := l90 ' 6 V <B := 99,214 

co B := 0.008 tcB :=0J288 a^ :=4.S98 

a m := 1.245- 10~ 2 , a 2B -=2-8* l0 ~ 6 a 3B : = ~ 2 - 703 ' 10 ~ 9 

C ^ :=4.184. (a (B + a m .T+ ajB'-T 1 ^ a 3B -*T 3 ) - R ^ 

It ie expected that similar coefficients will be programmed for a number of 
other inorganic gases. 
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tabulate the coefficients of equation 2. 



It is convenient to „„ 

f0 := 0.1445 I / co : = 6 - 3710 • • I- 



& 1 



:= 0.331 



f 2 := -0.1385 E2 



5 f 4 :=-6.07-l0 
2. Mixture properties 



-4 



£4 



:=-8- io 



, i_ m -fK known amounts of nitrogen and 



10 



15 



20 



25 



T cBH : ~ y T cfi" T cK 



3! 



z cAH := ^'^cA +z cH) 



v d3H ""f 



*cAB := '|' (r c^ +t cB) 
3. Yirial Coefficients and derivatives 
Let 



<B 



B H 



,V<H 



30 



m hvdrocarbon mixture are given by 
The virial coefficients for the unknown hydrocar 
(units of cc/mole) 2 1 

\3 f T <«\ 



cH 
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B2 H 



The virial coefficients for pure nitrogen are 

(f 0 +fiO-A) + fl-(^) + < f 2 +g 2-A)-(-^) - 

r /T/-a.\ 2 



V cA 
B A • 



(f 3 



B2 A :=— i2.f 1 -fe)+6.Cf2+g 2 ^A)'(^) - 

+ 12.(f 3 + g 3 - A)-(-^J + A) " V~/ 

The virial coefficients for pure methane are 



B1 B * 



B2 B 
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VcAH 

Bl xjl := 

z cAH 



. V cAH 
B2 AH •= — 



virial coefficients for the nitrogen-hydrocarbon mix are 

+ >3( f 3 +. g 3.0 ah)' (c-f 5 ) + -*< f 4 + * ^.AH)- ^—J 

2- f f t-|r-J + 6 ( f 2 + g 2 « AH)' \~~t~7 - 

' /X \ 3 /^cAHj 
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+ I2(f 3 +g 3 -a ah)^-tT-; +«V4 t M-a«^ t y 
The virial coefficients for the nitrogen-methane mix are 

B M :=^- f < f 0 + « <f« AB ) + f I' fer) + < f 2 + « I'* « );(^) ~ 

ICAB /T I 3 f T cAB\ 8 

+ (f 3 +g 3 o>AB)^- T ^j +a4+84'«^rr"j 



V cAB 
*cAB 



B2 ar 



The overall 



<f 0 +g 0 -<a ab ) 1 1" Y*~ t" . 
+ <f 3 + g 3 -*AB>-(^) +( f 4+£4' 

- f r(^) *« ab);(^) - . 

+ -3.<f3 + E 3 .a AB)'^) +-8 ! Cf 4 + £4-« AB)-^] 

- 2.f r (^^)+6(f 2 +g 2 -«AB)-(-^) 

/T \ 3 ^ T oAbV 
+ l2-(f 3 + E 3 -« ABV (-^J + **< f * + ^ 4^ AB)' 



+ l2-(f 3 +g 3 ' 

virial coefficients for the unknovm hydrocarbon 
Bx := B h-X H 2 + B a-X a 2 + 2 ' B AH" X lT X A 
Brf := Bl tf X + Bi^X A a * AlT X ir X - A 
Bx2 := B2 jfX H 2 + B2 a'X a 1 + 2 ' B2 AIT X lt X A 



are 
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and the overall virial coefficients for the known mixture are . 

Bl := Bl B -X B 2 + Bl A -X A 2 + 2 : Bl AB ?^ A -X B 
B2 := B2 B - X B 2 >HB2 iA - X A 2 + 2- B2 X A - X B 

The overall molecular weights, unknowa and known, are ! . 

mx := mH X jjH-m A X A mf:=.m, B -X B + m A -X A 

The gammas are \ f 

C vx— r X H + r X A T ** C w v 

^ gas gas ' , ,■ ■ 

CyB <=vA C v +l 

c Y :=_--x B + - — x A ' r >=-c7" 

gas gas v 

The molar volumes are 

R E as' T ■ .. R gas' T +p 

VX " = O.t01325-P +BX V '"0.101325-P 

where 0.101325 is the conversion of cc*atm to Pa*m A 3 

Finally, cx is the unknown hydrocarhon sound speed, which is a function of 
mH, and c is the sound speed of the known mixture. 



cx 



c 



cx = 423.0542 c t= 423.0542 



We choose the value mKE .015873 14 to make, c=cx. The unknown 
molecular weight is mx = 0.0176941 and the actual molecular weight is 

m « 0.0178385 for an error of <= '* TIie compressibility is 

m 

t :=ff Bx ' ai0132S,P wxthavalueof r«=0^704 
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The foregoing equations thus allow the determination vof sound speed , from 
molecular weight, and provide all ine : necessary relationships for eliminating or 
determining any unknowns if me temperature and pressure are known. Applicant applies 

5 this property of the molecular weight to determine the actual molecular weight by simply 
iteratively setting a test molecular weight, calculating the predicted sound speed, 
comparing to the actual sound speed, and using the error to set a new predicted 
molecular weight. The calculations are automated, and the sequence of test molecular 
weights quickly converges to the actual average molecular weightof the hydrocarbon 

10 portion of the measured gas. Using stored tables of the relationship between molecular 
weight and each of the gas properties, a computer program is readily written in BASIC 
or other suitable programming language to implement an iterative molecular weight 
determination as described above. 

15 i n a prototype operating embodiment, the processor is adapted to receive input 

data specifying the temperature, pressure and percentage of each of the above four 
standard inorganic gases, and to then calculate the speed of sound resulting from a 
mixture of these inorganic gases with a hydrocarbon mixture having an unknown average 
molecular weight. The processor follows the iteration procedure described above to 

20 estimate the actual average molecular weight using the foregoing empirical relationships 
and mixing rules under programmed control. This iteration procedure yields a value of 
the molecular weight which, if assumed for the hydrocarbon portion of the fluid in the 
conduit, accurately predicts the measured speed of sound. Discussion of representative 
calculations follows. 



25 



DISCUSSION AND EXAMPLES 
The molecular weight of an unknown combination of hydrocarbon gases in 
calculated from the speed of sound in the gas mixture. Applicant found this to be 
possible for a homologous series, such as the hydrocarbon gases, because the gas 
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properties which influence the speed of sound are also approximate functions of 
molecular weight. The .virial equation of states used to calculate the speed of sound in a 
gas under. non^ideal conditions, n For intermediate pressures (up to 20^50 bars), only die 
second virial coefficient is needed. This , is estimated from the gas critical data. In 
addition, data on molecular ; weight, :;W ^ and pressure are used. 

The virial equation for sound speed in a pure substance is illustrated by 
comparing the calculated (solid; lines) and measured (open circles and diamonds) speed of 
sound ethylene- This data is summarized. in Table I, and is graphically plotted in 
Figure 8. '■'•<..-: 



TABLE I 



15 



20 



25 



Sound Speed in Ethylene 
-20°C calculated . measured 
pressure (virial eqn.) <D 
(bar) m/£ in/s — 



0.00 
1-89 

9.01 

12.46 

16.63 

20.60 

23.65 

25J22 



309.60 
306.00 
300.08 
29L23 
282.94 
27157 
259.40 
249.10 
242.72 



308.46 
305.06 
299.42 
29L29 
283.68 
273.45 
26236 
253.08 
246.76 



100°C 
pressure 
(bar) 



calculated 
(virial eqa.) 
m/s 



measured 
(1) 
m/s 



0.00 

6.73 

10.74 

2SJ24 

28J54 

S7J94 



363.16 
359.42 
356.02 
346.94 
S42JBQ 
335.00 



363.77 
360.21 
357.08 
349.42 
346.28 
341.01 



The virial equation for the speed of sound inmixtures of known composition is 
illustrated by comparing the calculated (solid lines) and measured (open circles and 
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diamonds) speed of sound m two different nawrd ; gases, shown in Chart II and Chart III 
of Figures 9A and 9B, respectively. The corresponding data points are tabulated in 
TABLE H and TABLE III, respectively. As seen in tiiese examples, the measured speed 
of sound begins to deviate from the calculated speed of sound ^ pressure increases above 
about 20 bar. ^: >, v .. ; 



TABLE II TABLE ffl 





Press. 


calc. 


meafe. 






Press.. 


. calc. 


me as. 


Natural Gas 1 


(bar) 


(virial) 


(2) 




Natural Gas 2 


(bar) 


(virial) 


(2) 


composition 


0°C 


m/s 


m/s 




composition 


0«C 


cn/s 


m/s 


CI 81-20% 


0.O1 


405.14 


404.10 




" CI 84.60% 


0.01 


408.73 


407.60 


C2 1.05% 


6 


403.41 


402.40 




C2 1.46% 


6 


406.81 


405.40 


C3 0.04% 


10 


401.62 


400.90 




C3 3.47% 


10 


404-84 


403.40 


C4+ 0.04% 


20 


397.90 


398.10 




C4+ 0.06% 


20 


400.72 


400.10 


C02 0.4$% 


SO 


393.99 


895.80 




C02 0.69% 


30 


396.86 


397.40 


N2 17.17% 


40 


389.85 


894.10 




N2 13.84% 


40 


39L76 


395.30 


SO°C 










S0°C 








0.01 


423.98 


423.70 






0.01 


427.63 


427.40 




6 


422^5 


42250 






6 


426.46 


426.00 




10 


42190 


422.10 






10 


425.26 


424J90 




20 


419.76 


420.70 






20 


422£0 


423.00 




SO 


417.63 


419.70 






30 


420X5 


42L70 




40 


416X3 


419.10 






40 


417.61 


420.80 



20 



As noted above, using the relationship between molecular weight and each of the 
critical properties, and heat capacity, applicant is able to calculate the average molecular 
25 weight when only the temperature, pressure and speed of sound are measured. 

The method may be summarized as follows. All the hydrocarbon gases are 
considered to be a single component with an (arbitrarily assigned) molecular weight. A 
sound speed measurement is made. By using the virial equation and the correlations of 
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gas properties with molecular weight, a theoretical value for the speed of sound is also 
calculated. J The calciilaied value is compared to theobserved speed of sound,. and if 
differcn t, a new m^iecular weight value is chosen, . and, the calculations repeated, Thjs 
process is iterated untiTa niolecuiar weight is found which gives a calculated speed of 
sound equal to the observed speed of sound. Known amounts of inorganic components 
are considered as additional components of the mixture and are accounted for by using 
the mixing rules as above for virial coefficients. 

This virial method covers a temperature range of about -50°C to 300°C and a 
molecular weight range of about 2 to 130 g/ mole. A pressure range cannot be explicitly 
given, and will depend to some degree upon the particular gases involved. The method 
appears to produce satisfactory results at pressures up to 20 bars in all cases tested so for. 

RESULTS FOR MOLECULAR WEIGHT CALCULATIONS 

Errors resulting from application of the virial method to number of hydrocarbons 
at one atmosphere and 100°F were tabulated. Figure 10 compares the errors of this; 
method with the errors resulting from a currently used empirical equation. As shown, 
the present method yields more accurate results. As seen in Figure 10, the average error 
in calculated molecular weights for this set of forty -seven hydrocarbons plus hydrogen 
was 0.9%. 

A further test using experimental sound speed data was made on four natural 
gases. Here, the virial method was used to calculate a molecular weight from the 
reported sound speed/and this was compared to the molecular weight determined from 
the reported composition: For the purpose of calculation, the amounts of inorganic gases 
in each natural gas were taken as known. These results are shown in Table IV. 
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TABLE IV 



Natural Gas 1 MWc 
composition 

CH4 8L20% 


18.41= 0°C , 30*0 
Press. calc. calc. 
(bar) MW error * M W error 


C2H6 1.05% 
C3H8 0.04% 
C4H10* 0.04% 
O02 0.49% , 
N2 17.17% 1 


: 0.01 18.36 -4).27% , 1&32 -0.49% 
6 18.35 -0.33% 18130 -0.60% 
10 18J32 -0.49% 1&28 -0.71% 
20 18^5 -0.87% 1&22 -1.03% 
30 18.14 -1.47% 18.14 -1.47% 






Natural Gas 2 MW<= 
composition 

CH4 84.60% 
C2H6 1.46% 
C3H8 0.06% 
C4H10* 0.06% 
C02 0.59% 
N2 13.84% 


18.05 cpc sore 
Press. calc. calc. 
(bar) MW error MW. error 


0.01 18.00 -0.28% 1755 -0.65% 
6 18.01 -0.22% 1757 -0.44% 

,10 18.01 -0.22% 1756 -0.60% 
20 1754 -0.61% 1751 -0.78% 
30 17.82 -1.27% 17.83 -1.22% 






Natural Gas S MW*= 
composition 

GH4 82,62% 
C2H6 3.47% 
C3H8 0.76% 
C4H10* 0.39% 
O02 1.09% 
- N2 11.67% 


18.61 OXJ ; S0°C 
Press. . calc. calc. 
(bar) MW error MW error 


0.01 18.64 0.16% 18j68 -0.16% 
6 18.63 0.11% 18^7 -0.21% 
10 18.60 -0.05% 18-65 -0.82% 
20 lfl-60 -0.69% 18-48 -0.70% 
SO 18-36 -1.34% 18.38 -1.24% 






Natural Gas 11 MW«= 
composition 

CH4 83.44% 
G2H6 9.19% 
C3H8 3.67% 
C4H10* L18% 
C02 1.91% 
N2 0.71% 


19-45 0X3 SO-C 

Press. calc. calc. 
: (bar) MW error- MW error 


0.01 19.74 1.49% 19.71 1.34% 
6 19.72 1.39% 19168 L18% 
10 19.68 1.18% 19.64 0.98% 
20 19.64 0.46% 19.64 0.46% 
SO 19.33 -0.62% 19-39 -0.31% 



The average error shown in TABLE IV for the molecular weights of these natural 
gases was 0.8%. 
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The virial method described here uses nonlinear correlations of gas properties 
with molecular weight. Because of this nonlinearity , a wide distribution of molecular 
weights in a mixture will give an error in the properties predicted from the average 
molecular weight. In particular, an increased error occurs when a mixture contains 
hydrogen. In a Monte Carlo test the average errpr for mixtures containing hydrogen and 
four other randomly selected hydrocarbons was 2.3%. The worst error that could be 
produced by hydrogen m V a mixture was found 6.3%, illustratively for a mixture of 85% 
hydrogen and 15 % hexane. Distributions of this range are not commonly encountered, 
and this d^trib^^ of the calculate will not affect the overall accuracy or 

utility of the method for practical uses. 



RESULTS FOR COMPRESSIBILITY CALCULATIONS 

A significant advantage of applicant's virial method is that it also yields the 
compressibility factor of the gas mixture. This parameter is important for calculating 
n^s flow, particularly at pressures over one atmosphere. For example, for methane at 
20 atmospheres, using the compressibility factor increases the calculated mass by 3.4% . 
For heavier hydrocarbons, the increase would be even larger. Compressibilities z 
calculated for six gases (solid lines) are shown below in Figure 12 and are compared with 
handbook values (dark circles). 

Applicant's method actually gives the slope of compressibility with pressure, i.e., 
a straight line on a graph of compressibility vs. pressure. This works well in the 
intermediate pressure range for most gases, as shown above. While compressibility of 
nitrogen is nonlinear in this range, and is not predicted well, advantageously the overall 
at 40 bars is not increased if the calculated nitrogen compressibility factor were 



error 
used. 
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Thus, the invention accurately determines the average molecular weight of 
component gases of interest present in a gas, using only parameters such as temperature, 
pressure and sound speed that are directly and readily measured by ultrasonic 
interrogation or using simple sensors. The molecular weight may. then be used to 
identify upstream leaks, or sources or events giving rise to the detected hydrocarbon 
components; or may be used to calculate gas caloric value for furnace or energy capture 
feed streams; or to otherwise monitor or control flue or feed stream gases of otherwise 
unknown or poorly controlled composition. ; : : , , 



10 The invention being thus disclosed, further variations and modifications will 

occur to those skilled in the art, together with adaptations to known process control, 
monitoring or detection systems. All such variations, modifications and adaptations are 
considered to be within the scope of the invention, as set forth herein, and defined by the 
claims appended hereto and equivalents thereof. 



What is claimed is: 



20 
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Claims :" ' r * 

1. An ultrasonic measurement system tor aeiemumug a yi 

a conduit or WlfteVsucn^ r " • -r 

a plurality or tracer elements tanged for performmg. a signal path. : 
measurement of ultrasonic sighals mrough material, and 

a processor for receiving and processing such signals, wherein the processor is 
configured to determined average molecular weight of an unknown hydrocarbon mixture 
in said material. 

2 An ultrasonic measurement system according to claim 1. wherein the processor 
include a plurality of stored tables of critical constants of hydrocarbon mixtures as 
function of molecular weight of the mixtures. 

3. An ultrasonic measurement system according to claim 2, wherein the processor is 

configured to iteratively ; 

i) set a hypothetical hydrocarbon component molecular weight 

ii) look up the critical constants of the hypothetical molecular weight, 

iii) apply virial equations and mixing rules to determine a predicted sound 
speed for said hypothetical molecular weight hydrocarbon component together with one 
or more known inorganic components of said fluid material, and 

iv) repeat said steps ii-iii until the hypothetical molecular weight predicted 
sound speed is substantially equal to the measured sound speed thereby determining the 
average molecular weight of said unknown hydrocarbon mixture component. 

4 An ultrasonic measurement system according to claim 3, wherein said processor 
is configured for user input of known quantities of one or more gases selected from 
among nitrogen, hydrogen, carbon dioxide, sulfur dioxide and other inorganic gases. 
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5. A system for measurement or control, such system comprising 

a sound speed detector for measuring: sound speed in a fluid, the fluid including a 
background gas and an unknown fixture of one |r more hydrocarbons having an average 

molecular weight r jr > > 

5 a processor operative with the measured sound speed to determine the average 

molecular weight by iterative calculation with virial coefficients and mixing rules for the 
unknown mixture and background gas 

6. A system for measurement or control, such system comprising^ 

10 a sound speed detector for measuring sound speed in a fluid, the fluid including a 

background gas and an unknown mixture of one or more hydrocarbons having an average 
molecular weight 

a processor having a stored table of equivalences relating average molecular 
weight of a mixture of gases with virial coefficients, and operative with a measured 
15 sound speed to determine the average molecular weight by iterative calculation and 
thereby determine at least one critical property of the fluid. 
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Sound Speed in Nat. Gas #1 Sound Speed in Nat Gas 




pressure (bar) pressure (bar) 

Figure 9A Figure 9B 
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Oalc. and Meas. z in Acetylene 
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Figure 12(c) 

Calc. and Meas. z in Nitrogen 
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Gale and Meas. z in Carbon Dioxide 
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Calc. and Meas. z in Hydrogen 
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